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Abstract 

Acetylcholinesterase resistance has been well documented in many insects, including several mosquito species. We tested 
the resistance of five wild, Chinese strains of the mosquito Culex pipiens quinquefasciatus to two l<inds of pesticides, 
dichlorvos and propoxur. An acetylcholinesterase gene (acel) was cloned and sequenced from a pooled sample of 
mosquitoes from these five strains and the amino acids of five positions were found to vary (VISSM, G247S, A328S, A391T, 
and T682A). Analysis of the correlation between mutation frequencies and resistance levels (LC50) suggests that two point 
mutations, G247S (r^ = 0.732, P = 0.065) and A328S {r^ = 0.891, P = 0.016), are associated with resistance to propoxur but not 
to dichlorvos. Although the V185IV1 mutation was not associated with either dichlorvos or propoxur resistance, its RS 
genotype frequency was correlated with propoxur resistance {r^ = 0.815, P = 0.036). And the HWE test showed the A328S 
mutation is linked with V185l\/1, also with G247S mutation. This suggested that these three mutations may contribute 
synergistically to propoxur resistance. The T682A mutation was negatively correlated with propoxur (r^ = 0.788, P = 0.045) 
resistance. Knowledge of these mutations may help design strategies for managing pesticide resistance in wild mosquito 
populations. 
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Introduction 

Acetylcholinesterase (AChE, EC 3.1.1.7) is a key enzyme in the 
nervous system of both vertebrates and invertebrates that 
terminates nerve impulses by catalyzing tlie hydrolysis of the 
neurotransmitter acetylcholine (ACh) released from the presynap- 
tic membrane [1]. The inhibition of AChE by organophosphate 
and carbamate insecticides leads to the desensitization of the ACh 
receptor, thereby blocking nerve signal transmission. Organo- 
phosphates and carbamates have structures analogous to ACh and 
inhibit AChE competitively at the active site. Hydrolysis of these 
pesticide compounds retards the reactivation of the enzyme or 
inactivates it [2]. The extensive use of organophosphate and 
carbamate insecticides has resulted in the development of high 
levels of resistance to them among insects [3,4,,5,6]. 

Acel is the key AChE gene in insects. Several studies have 
found evidence that a point mutation in the ace 1 gene is associated 
with resistance to organophosphate and carbamate pesticides. This 
point mutation changes the structure of AChE making it 
insensitive to these insecticides. The first report of this mutation 
conferring insecticide resistance was in the two-spotted spider mite 
in 1964 [7]. Subsequent studies have demonstrated that many 
insect species have developed resistance to organophosphate and 



carbamate pesticides through decreased sensitivity of AChE [8], 
including many mosquito species, such as Anopheles gambiae [9], Cx. 
pipiens [10,11], Cx. pipiens quinquefasciatus [12], Cx. titaeniorhynchus 
and Cx. vishmi [13]. However, so far, only three acel mutations, 
G119S, F331W and F290V [T. califomica numbering) 
[13,14,15,16], have been confirmed to be involved in such 
resistance in mosquito species. Determining the mutations that 
confer resistance to specific pesticides is important to designing 
effective strategies for managing pesticide resistance. Cx. pipiens 
quinquefasciatus is the main mosquito species in urban environments 
in southern China and one of the most studied in terms of 
insecticide resistance. We here report the results of an investigation 
of mutations in the acel gene in five wild Chinese populations of 
Cx. pipiens quinquefasciatus. Knowledge of these mutations may have 
practical benefits for reducing pesticide resistance in this species. 

Results 

Resistance of the Five Mosquito Populations to 
Dichlorvos and Propoxur 

LCso values of the five different populations ranged from 0.266 
to 1.67 ppm for dichlorvos, and from 0.279 to 1.27 ppm for 
propoxur (Table 1). The HC strain had the lowest LCso and was 
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the most susceptible to both dichlorvos and propoxur. The SF 
strain had an LC5() to dichlorvos of 1.67 ppm and was 17.6 times 
more resistant to dichlorvos than the laboratory strain (LC50 
0.095 ppm). The GN strain had an LC50 to propoxur of 
1.27 ppm and was 11.0 times more resistant to propoxur than 
the laboratory strain (LC50 0.1 15 ppm). The HP strain was 7.89 
times more resistant to dichlorvos, and 4.62 times more resistant 
to propoxur, than the laboratory strain. The QB strain was 13.1 
times more resistant to dichlorvos, and 5.20 times more resistant 
to propoxur than the laboratory strain. 

Identification of Acel Mutations 

To identify mutations in the ace 1 gene, the cDNA of a pooled 
sample of mosquitoes from each of the five populations was 
cloned and sequenced. Five mutations (V185M, G247S, A328S, 
A391T, and T682A) in the pooled acel gene were identified 
(Figure 1), and the sequence was deposited in GenBank under the 
accession number KF680946. Note that this identification of 5 
mutations does not imply all occur in the same acel gene. The 
V185M mutation was GTG to ATG, the G247S mutation was 
GGC to AGC, the A328S mutation was GCG to TCC, the 
A391T mutation was GGC to ACC, and the T682A mutation 
was ACA to GCA. 

Polymorphism of the Acel Gene in Natural Population 

1. Determination of the allele frequencies. The allele 
frequencies of each mutation were determined by specific PGR 
amplification using the primers Gx-ace2-F, Gx-ace2-R and Gx- 
ace3-F, Gx-ace3-R on the cDNA obtained from individual 
mosquitoes. Genotypes of each mosquito in each population was 
determined by sequencing, and mutation frequencies (R%) 
computed (Table 2). We can see from Table 2 that the 
V185M, A328S and T682A mutations were present at different 
frequencies in all five strains. However, the A391T mutation was 
only found in the HP and QB strains, and the G247S mutation 
was found in all but the HG strain. 

2. Hardy— Weinberg Equilibrium (HWE) test and genetic 
linkage analysis of the mutations. The results of GENE- 
POP software analysis of HWE and genetic linkage of the 
acetylcholinesterase gene mutations are shown in Tables 2 and 3. 
The HWE test indicates the QB and GN populations have a 
heterozygote deficit with respect to the T682A mutation (P< 
0.05), and the HP population a heterozygote excess with respect 
to the A391T mutation (P<0.05). Mutations in all other 
populations did not deviate from the HWE and none of the five 
mutations deviated from the HWE across all populations (P> 
0.05). 

Results of linkage disequilibrium analysis of the five mutations 
are shown in Table 3. Evidence of linkage disequilibrium was 
found for V185M with respect to the A328S and A391T 
mutations (P<0.05), The G247S and A328S mutations' hnkage 
disequilibrium P-value was 0.0821, only slightiy above 0.05. This 
suggests that these two mutations might exist in the same gene. 
Our sequencing data indicated that that these two mutations do 
indeed occur in the same acel gene in some mosquitoes. But the 
conclusion had to be confirmed by more data. All otiier gene 
polymorphism was randomly distributed. 

Correlation of Resistance with Mutation Frequencies 

The correlation between resistance to dichlorvos and propoxur 
and the frequencies of four mutations (V185M, G247S, A328S, 
T682A) are shown in Figure 2 and Table 4. The four mutations' 
frequencies were all not significantly correlated with dichlorvos 
resistance. Although the frequency of the V185M mutation was 
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uncorrelated with propoxur resistance (Figure 2 A), its RS 
genotype frequency was (r^ = 0.815, P = 0.036) (Figure 2 B). The 
correlation between the frequency of the G247S mutation and 
propoxur resistance was close to significance (r^ = 0.732, 
P = 0.065), and there was a significant linear relationship between 
the frequency of the A328S mutation and propoxur resistance 
(r^^ 0.891, P = 0.016) (Figure 2 C, D). The frequency of the 
T682A mutation was negatively correlated with propoxur 
(r^ = 0.788, P = 0.045) resistance (Figure 2 E). 

3D Models of Mutations and Structural changes at the 
Catalytic Site 

A 3D model was made of the Cx. pipiens quinquefasciatus ace 1 gene 
sequence allowing the location and structure of four mutations to 
be visualized (Figure 3). The V185M and A391T mutations are 
distant from the active site of the enzyme-catalytic triad (S327, 
H567, £453; S200, H440, £327 in T. califomica) (Figure 3A, B). 
The other two mutations, G247S and A328S, are close to the 
catalytic site (Figure 3C, D) and could therefore potentially affect 
the binding between ACh£ and its substrates (Ach: ZINC3079336 
and propoxur: ZINC 1590885). Figure 3£-H illustrates the change 
in amino acids and H-bonds associated with the G247S and 
A328S mutations. These two substitutions change the amino acids 
present at catalytic sites removing the two H-bonds (S327(8) 
0^,-03, S327(8)Oy-04) between ACh£ and Ach (Figure 3E, F) 
and reducing the three H-bonds between ACh£ and propoxur 
(G247(4)-013, S327(8)0^-011, H567(14)-NH27) to one 
(S327(10)Oy-NH27) (Figwe 3G, H). Hence, these two mutations 
could have a major effect on the catalytic activity of the AChE 
enzyme. 

Discussion 
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The indiscriminate use of insecticides over more than half a 
century has resulted in high levels of insecticide resistance in many 
mosquito species [13,17,18]. We tested the resistance of five 
Chinese Cx. pipiens quinquefasciatus populations to dichlorvos and 
propoxur. Our results show that, compared to a laboratory strain, 
these five populations displayed a 2.80- to 17.6-fold resistance to 
dichlorvos and 2.43- to 11.0-fold resistance to propoxur. The 
frequent use of these insecticides has created an intense selection 
pressure for traits that confer resistance to them, such as changes 
in behavior, epidermal structure, metabolic enzymes and target 
site mutations. Resistance may be conferred by any one, or more 
than one of these mechanisms. Osta et al (20 1 2) found that the 
dramatic reduction in the frequency of the G119S [T. califomica 
numbering) mutation in Culex pipiens mosquitoes was probably due 
to the increased use of pyrethroids over organosphosphate 
insecticides [19]. Therefore, alternating between different kinds 
of insecticides is one way of minimizing the development of 
resistance to any one kind. 

We used cloning and sequencing to identify five point mutations 
in the acel gene of Chinese Cx. pipiens quinquefasciatus. HWE tests 
suggest that these five mutations do not deviate from the HWE 
across all populations. However, the tests also indicated that the 
QB and GN populations were deficient in heterozygotes with 
respect to the T682A mutation and that HP population had an 
excess of heterozygotes with respect to the A391T mutation (P< 
0.05). Further work will be required to determine the reasons for 
these departures from the HWE. Linkage disequilibrium analysis 
indicated significant linkage between the V185M mutation and 
the A328S and A391T mutations. Although linkage between other 
mutations was statistically insignificant, that between the G247S 
and A328S mutations was nearly so (P= 0.0821). Our sequencing 
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Figure 1 . Alignment of nucleotide and amino acid sequences of 

Cx. pipiens quinquefasciatus. Cxql is the template nucleotide 
sequence (no amino acid mutation) and Cxq2 the mutant nucleotide 
sequence. Nucleotides are numbered on the first line, amino acids on 
the second. The five mutations are shown in the black frames. 
doi:10.1371/journal.pone.0095260.g001 
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Table 3. P-value for linkage disequilibrium of each pair of loci across all populations (Fisher's method). 



Locus pair 




df 


P-Value 


V185M & G247S 


11.237 


8 


0.1887 


V185M & A328S 


23.804 


10 


0.0081 


G247S & A328S 


13.988 


8 


0.0821 


V185M & A391T 


7.5840 


2 


0.0226 


G247S & A391T 


3.4992 


2 


0.1738 


A328S & A391T 


5.8691 


4 


0.2091 


V185M & T682A 


4.8208 


10 


0.9028 


G247S & T682A 


2.4273 


6 


0.8765 


A328S & T682A 


14.160 


10 


0.1658 


A391T & T682A 


1.1946 


4 


0.8790 



doi:l 0.1 371 /journal.pone.0095260.t003 



results suggest that these two mutations occur within the same 
acel gene in some mosquitoes but further work is required to 
confirm this hypotliesis. 

These results are the first report of the V185M mutation in Cx. 
pipiens quinquefasciatus. Although there was no apparent correlation 
between the frequency of this mutation and resistance to 
dichlorvos and propoxur, the frequency of its RS genotype was 
significantly correlated with propoxur resistance (r^ = 0.815, 
P = 0.036). Although the 3D model indicates that V185M is 
located far from the active site, the positive linear relationship 
between its RS genotype frequency and propoxur resistance, and 
its apparent linkage with the A328S mutation suggest that it may 
be involved in propoxur resistance. Of course, we cannot rule out 
the possibility that insecticide resistance involves multiple dupli- 



cation of the acel gene. Further research needs be required to 
determine how this might affect the catalytic center. 

Our results (Figure 2, Table 4) suggest that the G247S mutation 
is not associated with propoxur resistance and that the A328S 
mutation is. The G247S mutation corresponds to G119S in T. 
califomica which has been associated with insecticide resistance in 
mosquitoes by several authors [20,21]. The Gl 19 position is part 
of the oxyanion hole (Gl 18, Gl 19, and A20 1 in T. califomica), close 
to the catalytic Serine (S200) where a G to S substitution would 
reduce accessibility to inhibitors and substrate by steric hindrance. 
SI 19 is close enough to the catalytic residues to alter the 
presentation of inhibitors and substrates. This could be the reason 
this mutation confers resistance to some insecticides [22,23]. 
Although the correlation between the frequency of the G247S 
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Figure 2. Linear regression of the relationship between resistance levels (LC50) and mutation ratios. Resistance levels to propoxur are 
plotted against the ratios of V185IV1 (A), the RS ratio of V185M (B), G247S (C), A328S (D), and T682A (E). 
doi:1 0.1 371 /journal.pone.0095260.g002 
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Table 4. The analysis results of correlation between propoxur LC50 and mutation frequencies. 





Mutations 


Insecticide 


R (95% CL^ 




P 


Significance^ 


V185M 


Propoxur 


0.647{-0.549,0.974) 


0.419 


0.238 


No 


V185M (RS%) 


Propoxur 


0.903(0.101,0.994) 


0.815 


0.036 


Yes 


G247S 


Propoxur 


0.855{-0.1 10,0.990) 


0.732 


0.065 


No 


A328S 


Propoxur 


0.944(0.366,0.996) 


0.891 


0.016 


Yes 


T682A 


Propoxur 


-0.887(-0.993, -0.023) 


0.788 


0.045 


Yes 



^CL = confidence limits. 
^a = 0.05. 

doi:l 0.1 371 /journal.pone.0095260.t004 



mutation and propoxur resistance was not statistically significant 
(P = 0.065), numerous prior publications have reported such an 
association [9,10,12] and noted that this mutation is often 
combined with other mutations in resistant strains. Therefore, 
we suspect that G247S probably is involved in propoxur 
resistance. We may have failed to detect a significant correlation 
between the frequency of this mutation and resistance because of 
its low frequency in our sample, which could be because most 
mosquitoes carrying it were heterozygotes. Furthermore, the 
resistance conferred by this mutation may be nearly recessive 
under certain bioassay conditions [24]. 

The A328 position corresponds to the A201 position in T. 
califomka, which is located within the active gorge of the enzyme, 
close to the catalytic site,and is a part of the oxyanion hole. Li et al 

(2009) also found the A328S mutation in Cx. pipiens pallms and 
made a three-dimensional model of AChE to visualize this 
mutation. However, they did not demonstrate a relationship 
between the A328S mutation and resistance [25]. Khajehali et al 

(2010) found the corresponding A201S mutation in Tetranjchus 
urticae Koch, and demonstrated that this was possibly involved in 
resistance to organophosphorus and carbamate insecticides [26]. 
Our results suggest that this mutation is involved in propoxur 
resistance (r^ = 0.891, P = 0.016). The linkage disequilibrium and 
sequencing results indicate that A328S and G247S mutations exist 
in a same acel gene in some mosquitoes, which suggests that they 
may work synergistically. The G119 and A201 positions {T. 
califomica numbering) are both part of the oxyanion hole, and 
could therefore both contribute an amide nitrogen to form bonds 
that could stabihze the enzyme-substrate complex. The substitu- 
tion of serine for glycine and alanine may change the conjunctions, 
conferring resistance to some insecticides [23,27]. We can see from 
Figure 3 that these substitutions could decrease the numbers of H- 
bonds between enzyme and substrate. H-bonds are the strongest 
force between molecules so a reduction in these could reduce 
enzyme-substrate stability and interfere with the catalytic reaction. 

The Gl 19S mutation was the first mutation found in mosquito 
vectors [10]. Previous studies indicated that this mutation would 
incur a high fitness cost [28], however, although the cost of 
resistance is often high at the beginning of selection when 
resistance is unstable, the cost reduces and resistance stabilizes 
with increasing duration of exposure to insecticides [29]. Other 
mutations can play an important role in this process. Mutero et al 
(1994) found that high levels of resistance were obtained by the 
combination of several point mutations [1] and Menozzi et al 
(2004) demonstrated that combining mutations could increase 
insecticide resistance in Drosophila melanogaster [30]. Our results 
show that the A328S mutation (A201S in T. califomica) may work 
synergistically with the Gl 19S mutation in the oxyanion hole. It's 
possible that the A328S mutation compensates for some of the 



fitness costs incurred by the G 1 1 9S mutation. This is a fascinating 
question but further in vitro assays are required to confirm this 
hypothesis. 

The A391T mutation was found only in the HP and QB 
populations, which had moderate LC50 values and in which it had 
a frequency of around 0.500. The genetic linkage analysis 
indicates a linkage between this mutation and V185M, however, 
in view of the small sample size further work is required to confirm 
this. The three-dimensional model revealed that the A391 
mutation is distant from the active site. This indicates that this 
mutation is unlikely to alfect catalytic activity and is probably not 
involved in dichlorvos and propoxur resistance. How this mutation 
developed and its function, if any, in pesticide resistance requires 
further investigation. 

Our results provide the first evidence of the T682A mutation in 
Cx. pipiens quinquefasciatus. The frequency of this mutation was 
negatively correlated with propoxur resistance (r^ = 0.788, 
P= 0.045). Fournier et al (1988) found that AChE in Drosophila 
melanogaster was composed of two, non-covalendy associated, 
polypeptides of 55 and 16 kDa. AChE is an amphiphilic protein 
linked to the membrane of neuronal cholinergic synapses via a 
glycolipid anchor at the C-terminal end of the 55 kDa polypeptide 
[31]. Nabeshima et al (2004) found an I697M replacement near 
the C-terminus (Ile701) in Culex tritaeniorhynchus, but considered 
that this was unlikely to be the cause of AChE insensitivity [32]. 
Our results also indicate that the T682A mutation is near the C- 
terminus of AChE, and that the frequency of this mutation is 
negatively correlated with propoxur resistance. Despite its negative 
correlation with resistance, it's possible that this mutation may 
change the C-terminus structure of AChE thereby reducing its 
attachment to the membrane and the stability of enzyme. We 
don't know whether this mutation works in combination with the 
other four mutations or not, or if its apparent negative relationship 
with resistance is related to fitness costs. 

In conclusion, we found five ace 1 gene mutations in Cx. pipiens 
quinquefasciatus that are correlated with propoxur, but not 
dichlorvos resistance. The V185M mutation was first confirmed 
in Cx. pipiens quinquefasciatus and may be involed in propoxur 
resistance. The allele frequencies of the G247S and A328S 
mutations were positively correlated with resistance. So the G247S 
and A328S mutations are also likely to confer propoxur resistance. 
The A391T mutation appears unrelated to dichlorvos and 
propoxur resistance and the T682A mutation appears negatively 
correlated with resistance to propoxur. Ideiitilyiiig the mutations 
that confer resistance to specific insecticides can inform the choice 
of insecticides for a given insect population, thereby reducing the 
development of resistance and improving the efficacy of control. 
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Figure 3. Three dimensional model of the ACHE of Or. pipiens quinquefasciatus based on the structure of T. californica (PDB: 3ZV7). The 

four mutations are shown as red, green and blue van der Waals spheres and the catalytic triad (S327, H567, E453; S200, H440, E327 in T. californica) is 
shown in yellow. A-D illustrates the four mutations. A shows the VI 85 and A391 positions and B the IVI185 and T391 mutations. C and D show the 
G247S and A328S mutations, and the catalytic triad. E-H shows changes in the enzyme-substrate complex; Ach (E, F) and propoxur (G, H) are shown 
in green and the H-bond as yellow dotted lines. Amino acids are marked with numbers. E1-14 (wild-type enzyme) are W212, G245, G246, G247, Y249, 
Y258, E326, S327, F416, Y456, F457, H567, G568, 1571 respectively; The two H-bonds were composed of S327 Ogamma and 03, S327 Ogamma and 
04. F1-14 (G247S/A328S mutant) are W212, F244, G245, G246, S250, G251, T252, L255, Y258, S327, Y456, H567, G568, 1571 respectively; G1-16 (wild- 
type enzyme) are W212, G245, G246, G247, Y249, Y258, E326, S327, W360, F416, Y456, F457, F527, H567, G568, 1571 respectively; The three H-bonds 
were composed of G247(4) NH and 013, S327(8) Ogamma and 01 1, H567(14) and NH27. H1-17 (G247S/A328S mutant) are W212, F244, G245, G246, 
S247, Y249, L255, Y258, E326, S327, W360, F41 6, Y456, F457, F527, H567, 1571 respectively. The only H-bond was composed of 5327(1 0) Ogamma and 
NH27. 

doi:1 0.1 371 /journal.pone.0095260.g003 
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2109bp 

5' 3' 

420bp 795bp 

* ► 

Cx-ace1-F Cx-ace1-R 932bp Cx-ace3-F Cx-ace3-R 

► 

Cx-ace2-F Cx-ace2-R 



Figure 4. Schematic diagram of the amplification of the acel gene. The complete sequence was 2109 bp (blacl<), and the three sections are 

indicated by red, blue and green arrows. 

doi:10.1371/journal.pone.0095260.g004 



Materials and Methods 

Statement of Ethical Approval 

No ethical approval was required as no regulated animals were 
used in this study. Pre-permission (May-September 2012) was 
granted for observation, collection and field research on mosqui- 
toes in Guangdong and Hainan Provinces, which was conducted 
as part of the Infective Diseases Prevention and Cure Project of the 
China National Ministry of Pubhc Health (No.2008ZX10004 and 
No.2012ZX10004219). All field studies on Cx. pipiens quinquefascia- 
tus were authorized by Guangdong and Hainan Provincial CDC 
Committees for Animal Welfare and Animal Ethics (address: 1 76 
Xingang West Road, Guangzhou, Guangdong province, and 44 
Haifu Road, Haikou, Hainan province, P. R. China). 

Mosquito Strains 

Specimens of Cx. pipiens quinquefasciatus were collected from five 
difierent field sites; Guangzhou Nansha (El 13°29'29.35", 
N22°48'4.13") and Haikou Poxiang (El 10°19'33.79", 
N19°59'55.07") in May 2012, and Haikou Changliu 
(El 10° 11 '50.36", N20°0'50.25"), Qionghai Boao 
(E110°34'57.13", N19°09'42.07") and Sanya Fenghuang 
(E109°26'54.38", N18°18'2.91") in September 2012. The suscep- 
tible strain had been reared in an insectarium for more than 10 
years without exposure to any insecticides. 

Cx. pipiens quinquefasciatus larvae were collected at each field site 
and reared to adulthood. Some wild caught female adults were 
frozen in liquid nitrogen for subsequent testing. 

Bioassay 

Bioassays were conducted by putting thirty late 3"' or early 4* 
instar larvae into pans containing 200 ml water. Measured 
quantities of insecticides were added to each pan using an 
automatic pipette according to the methods specified by the WHO 
[33]. Larval mortality was recorded 24 h after each treatment. No 
food was ofiered to larvae during bioassays. Larvae were 



maintained in the laboratory under a 14L:10D photoperiod, 
75% relative humidity and temperature of 26±1°C during 
bioassays. Bioassays of each insecticide were repeated three times. 
Statistical analyses were performed using SPSS software version 
13. 

Extraction of RNA and cDNA Synthesis 

Total RNA was extracted from specimens from each population 
with Trizol reagent (GBT) following the manufacturer's protocol 
and cDNA synthesized from the total RNA with cDNA synthesis 
kit (TaKaRa). The cDNA was stored at -20°C. 

PCR Amplification 

Gene specific primers based on the pubhshed insecticide 
resistant sequence of the Cx. pipiens quinquefasciatus acel gene 
(GenBank Accession No.:CQ,753634.1, this includes a G119S 
mutation related to propoxur resistance) were designed in NCBL 
Primer-BLAST and used to amplify the acel gene of each 
population. The acel gene is 2109 bp and is divided into three 
sections (Figure 4). The primers used are shown in Table 5. 

Cloning and Sequencing of PCR Products 

PCR products were purified using a universal DNA purification 
kit (TL\NGEN) and the purified products were ligated into the 
pEASY-Tl vector (TRANSGEN). The recombinant plasmids 
were cloned into Trans 1-Tl competent cells (TRANSGEN). The 
microbials were spread on LB solid medium (including ampicillin, 
X-gal, IPTG) and cultured overnight. White clones were selected, 
placed in LB liquid medium and cultured to turbidity. Positive 
clones were identified by PCR using Ml 3 forward and reverse 
primers and sequenced by Sangon Biotech [25]. Based on the 
discovery of clones, the genotype of individual mosquitoes was 
determined for each amino acid position by specific PCR 
amplification and sequencing. Li this procedure, a single 
mosquito's RNA was extracted and reversed transcribed to cDNA, 



Table 5. The primers used to amplify the Cx. pipiens quinquefasciatus acel gene. 



Primers 5'^3'Sequence length (bp) PCR parameters 

Cx-acel-F ATGGAGATCCGAGGCCTAAT 420 94°C,5 min; 94 C,30 s; 62=C,30 s; 72°C,1 min,35 cycles; 72=C,7 min. 

Cx-acel-R GCCCTTGTCCGTCGTTATG 

Cx-ace2-F CGGACCCACTGGTCATAACG 932 94°C,5 min; 94 C,30 s; 65=C,30 s; 72°C,1 min,35 cycles; 72=C,7 min. 

Cx-ace2-R ACCCTCCTCGGTGTTGCTG 

Cx-ace3-F CGCTTCAAGAAAACGGA 795 94°C,5 min; 94 C,30 s; 55=C,30 s; 72°C,1 min,35 cycles; 72=C,7 min. 

Cx-ace3-R TTAAATCTTGAACCGCGT 



doi:l 0.1 371 /journal.pone.0095260.t005 
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then amplified by specific PGR before sequencing. Calculated 
mutation fi-equencies were based on the sequencing results. 

Hardy-Weinberg Equilibrium (HWE) Test and Genetic 
Linkage Analysis of the Mutations 

The Hardy-Weinberg equilibrium (HWE) describes the theo- 
retical fi-equency of two alleles of a single locus in the absence of 
mutation and selection after one generation of random mating in 
an indefinitely large population ^\ ith dis('rc't(; g(;n("rations [34] . We 
used GENEPOP software to analyze the HWE and genetic linkage 
of mutations. 

Correlation of Pesticide Resistance with the Allele 
Frequency of Different Mutations in the Five Mosquito 
Populations 

The resistance (LG50) of the five populations to propoxur and 
dichlor\'OS was determined by bioassay and the allele frequencies 
of the various mutations were determined by gene specific 
amplification and sequencing as described above. The LG50 of a 
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